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ABSTRACT 

Polarization measurements provide direct insight into the nature of astrophysical processes. 
Unfortunately, only a few instruments are available for this kind of measurements at 7-ray ener- 
gies, and the sources need to be very bright. Gamma-Ray Bursts (GRBs) are ideal candidates 
due to their large flux over limited time intervals, maximizing the available signal-to-noise ra- 
tio. To date a few polarization measurements have been reported, claiming of a high degree of 
polarization in the prompt emission of GRBs but with low statistical evidence. 

We used the IBIS telescope on board the INTEGRAL satellite to measure the polarization 
of the prompt 7-ray emission of the long and bright GRB 041219A in the 200-800 keV energy 
band. We find a variable degree of polarization ranging from less than 4% over the first peak to 
43±25% for the whole second peak. Time resolved analysis of both peaks indicates a high degree 
of polarization, and the null average polarization in the first peak can be explained by the rapid 
variations observed in the polarization angle and degree. 

Our results are consistent with different models for the prompt emission of GRBs at these 
energies, but they favor synchrotron radiation from a relativistic outflow with a magnetic field 
which is coherent on an angular size comparable with the angular size of the emitting region 
(^1/r) . Indeed this model has the best capabilities to maintain a high polarization level, and 
to produce the observed variability. 

Subject headings: gamma rays: bursts — polarization — gamma rays: observations 



1. Introduction 

Gamma Ray Bursts (GRBs) are powerful 
flashes of 7-ray radiation appearing at random 
directions on the sky, and lasting from a fraction 
to hundreds of seconds. Some of them are firmly 
associated with supernovae of type Ib/c, and the 
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general picture is that bursts longer than 1 s are 
associated with the death of massive ( >, 30 M(7i) 
stars (for a recent review on GRBs, see Meszaroa 
20061 ). The isotropic equivalent energy radiated 



during the prompt phase is of the order of 10 to 
10^'^ergs and is believed to originate from a highly 
relativistic outflow (F > 100) ejected by the cen- 
tral source. The precise content of this jet, and 
especially its magnetization, as well as the details 
of the mechanism leading to the 7-ray emission 
are still uncertain. Models range from unmagne- 
tized flreballs where the observed emission could 
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be produced by relativistic electrons accelerated 
in interna l shocks propagating within the outflow 
(jRees fc M cszaros 1994), to pure electromagnetic 
outflows where the r adiated energy comes from 
magnetic dissipation ( Lvutikov 20061 ). Intermedi- 
ate cases with mildl y magnetized outfl ows are of 
course possible (e.g. ISpruit et al.ll200l[ ). Even in 
the case of an unmagnetized fireball, a local mag- 
netic field in the emission region, generated by 
the shocks, is necessary if the dominant process is 
synchrotron radiation from relativistic electrons. 
In the case of a mildly magnetized outfiow, such 
additional shock-generated magnetic field could 
also be present. The present situation is that nei- 
ther the global magnetization of the outflow, nor 
the local intensity of the field are well constrained 
by the observations. 

Polarization measurements of the prompt phase 
of GRBs can shed new light on the strength and 
scale of magnetic fields, as well as on the radia- 
tive mechanisms at work. Locally, the synchrotron 
emission has a high degree of linear polarization, 
that however could be reduced by relativistic ef- 
fects in a jet. In the case where the magnetic field 
is mainly transverse and highly ordered, i.e. has a 
coherence scale which is larger than the typical size 
~ R/T of the visible part of the emitting region, 
the detected signal can still be highly polarized. 
On the other hand, in the case of a random field 
or an ordered magnetic field parallel to the ex- 
pansion velocity, the polarization of the detected 
signal should vanish, except for the peculia r con- 
dition of a jet observed slightly off- axis (,Waxman 
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To 



date only a few polar ization measure- 

Coburn fc Boggd 



ments are available for GRBs. 



(|2003l ) reported a high degree of polarization, 
n = 80% ± 20%, for GRB 021206. However, 
successive reanalysis of the same data set could 
not confirm this claim, reporting a degree of po- 
larizatio n compatible wi th zero (iRutledge fc Foxl 
[2004- ,Wigger et all 12004 . IWillis et all (l2005l) re- 
ported a strong polarization signal in GRB 930131 
(n > 35%) and GRB 960924 (B > 50%), but 
this result could not be statistically constrained. 
GRB 0412 19 A is the longes t and brightes t GRB 



localized by INTEGRAL (IWinkler et all 120031) 



to da te (jVianello et al.l 120091 ). iMcGlvnn et al 



(|2007l ). using t he data of the INTE GRAL spec- 
trometer (SPI; IVedrenne et al. I l2003h . reported a 



high degree of polarization of the prompt emission 
(n = 68 ±29%) for the brightest part of this GRB. 

Here we report the polarization measurement 
of the prompt emission of GRB 041219A, per- 
forme d using the imager on board INTEGRAL^ 
IBIS dUbertini et al.ll2003l) . 



2. Observations and analysis 

IBIS is composed of two position sensitive 
d etector layers IS GRI (CdTe, 15-1000 keV; 
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Lebrun et al.l 120031). and PI CsIT (Csl, 200 keV- 
10 MeV: iLabanti et al.|[2003l). and can be used as 
a Compton polarimeter ( Lei et allll997l ). thanks 
to the polarization dependency of the differential 

cross section for Compton scattering 



2 sin^ 9 cos^ < 



(1) 

where Tq is the classical electron radius, Eq the 
energy of the incident photon, E' the energy of 
the scattered photon, 9 the scattering angle, and 
(f) the azimuthal angle relative to the polarization 
direction. Linearly polarized photons scatter pref- 
erentially perpendicularly to the incident polariza- 
tion vector. Hence by examining the scatter angle 
distribution of the detected photons 



A^(0) -5[l + aoCos2((/.-0o)], 



(2) 



one can derive the polarization angle, PA =^ 
00 — 7r/2 -|- nTT, and the polarization fraction H = 
Oq / aioo , where Oioo is the amplitude expected for a 
100% polariz ed source derived by Montecarlo sim- 
ulations (see Forot et al. 20081 ). 

IBIS is a coded mask telescope with a 29°x29° 
field of view at zero sensitivity, and a 8.5° x8. 5° 
central region where the sensitivity is maximal 
and uniform. GRB 041219A was detected at 
~3.15° off-axis. Here we consider the events which 
interacted once in the upper layer, ISGRI, and 
once in the lower layer, PICsIT, and whose recon- 
structed energy lies in the 200 keV-800 keV range. 
These events are automatically selected on board 
through a time coincidence algorithm, whose max- 
imal allowed time window was 3.8 /is during our 
observation. A light curve of the Compton events 
is reported in Fig. [1] 

As can be seen, the GRB is clearly detected, but 
due to its brightness, the available IBIS telemetry 
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Fig. 1.— Compton light curve of GRB 041219A. 
Bin size is 5 s. The two insets show a magnified 
view of the two peaks, binned at 1 s. The ana- 
lyzed intervals, are shown with dashed lines. P8 
is omitted for clarity. 

share was saturated, and hence not all the events 
could be sent to the ground. Images based on 
Compton events selected in the first and second 
peak show the bursting source respectively at a 32 
sigma level and 20 sigma level. An image made 
on the whole GRB is shown in Fig [21 and the 
position derived (aj2ooo = 00''24™24^ ^j20oo = 
+62°49'44" with an uncertainty of 1' at 90% c.l.) 
is consistent with the one derived from simultane- 
ous optical data ( Vestrand et al. 2005f) . 

To measure the polarizati on, we followed the 
same procedure described in Forot et al. ( 2008h 



that allowed to successfully detect a polarized sig- 
nal from the Crab nebula. To derive the source 
flux as a function of (p, the Compton photons were 
divided in 6 bins of 30° as a function of the az- 
imuthal scattering angle. To improve the signal- 
to-noise ratio in each bin, we took advantage of 
the TT-symmctry of the differential cross section 
(see Eq. [T]), i.e. the first bin contains the pho- 
tons with 0° < < 30° and 180° < (/) < 210°, 
etc. Then the chance coincidences (i.e. photons 
interacting in both detectors, but not related to a 
Compton event), have been subtracted from each 
de tector image follow ing the procedure described 



Forot et al.l (|2008l ). The derived detector im- 



ages were then deconvolved to obtain sky images, 
where the flux of the source in each bin is mea- 
sured by fitting the instrumental PSF to the source 




Fig. 2. — 200-800 keV Compton significance im- 
age of GRB 041219A. The source is detected at 37 
sigma level. The white cross indicates the positio n 



of the optical counterpart ( Vestrand et al. 2005h . 



and grey contours are traced at equal significance 
starting from 3 in steps of 6. 



peak. We finally fitted using a least squares tech- 
nique the polarigrams (see Fig. [3]) with Eq. [2] to 
derive oq and (po, and the errors on the parameters 
are dominated by the statistics of the data points. 
To evaluate the goodness of our fits, we computed 
the c hance probability (see Eq. 2 in iForot et al.l 



20081 ) that our polarigrams are due to an unpo- 
larized signal, and reported these values in Fig. 

m 

3. Results 

We analyzed the different portions of the GRB, 
focusing on the brightest parts. First we ana- 
lyzed the entire first and second peak, and then 
we performed a time resolved analysis: 36 inter- 
vals lasting 10 s, each one overlapping for 5 s with 
the previous one, have been analyzed over the 
whole duration of the GRB starting at 01:46:22 
U.T. until 01:49:22 U.T. The most significant non- 
overlapping ones (P6, P8, P28, P30) have been 
chosen for the polarization analysislj. The inte- 
gration times and the imaging significance of the 
chosen time intervals are reported in Table [1] and 



^Even if P9 is more significant than P8 (see Table [TJ, we 
prefer to discuss P8 because it is simultaneous to the SPI 
analysis. Anyway, P8 and P9 results are statistically con- 
sistent. 
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Fig. [TJ The azimuthal distributions of the GRB 
flux for the different time intervals are reported in 
Fig. [3 



First Peak P = O.S 



Second Peak P = ^^.IxlO" 




100 200 300 iOO 300 300 

Fig. 3. — Polarigrams of the different time inter- 
vals that have been analysed (see Table [ij . For 
comparison purposes, the curves have been nor- 
malized to their average flux level. The crosses 
represent the data points (replicated once for clar- 
ity) and the continuous line the fit done using Eq. 
m The chance probability of a non-polarized sig- 
nal is reported in each panel. 

As can be seen, no polarization signal could be 
found integrating over the whole first peak, and 
the upper limit is 4%. On the other hand, a modu- 
lated signal is seen in the second peak correspond- 
ing to n = 43 ± 25%. Integrating over smaller 
portions of the GRB, we measure highly polarized 
signals, especially in P8, P9 and P30. 

4. Discussion and conclusions 



Using SPI data lMcGlvnn et all (|2007l ) reported 
high polarization for GRB 041219A, but they re- 
stricted their analysis to the first peak of the GRB. 
They analyzed a 66 s interval within the first peak, 
and a 12 s interval on the brightest part of it start- 
ing at 01:46:54 U.T. The latter corresponds to our 
P8 interval and indeed our results are compatible 
at 1 sigma level, both for the polarization fraction 
and angle, for which they find n=68±29% and 
PA=70tloo (100 keV-1 MeV). On the other hand 
we cannot confirm their result on the broader 
66 s time interval (n=26±20% and ^^=70^2?=), 
starting at 01:46:21 U.T., which corresponds to 



Table 1: Polarization results for the different time 




intervals. 












Name 


start 




n 


PA 


Image 




U.T. 


U.T. 


% 


degrees 


SNR 


First Peak 


01:46:22 


01:47:40 


<4 




32.0 


Second Peak 


01:48:12 


01:48:52 


43±25 


38±16 


20.0 


P6 


01:46:47 


01:46:57 


22± 13 


121±17 


21.5 


P8 


01:46:57 


01:27:07 


65±26 


88±12 


15.9 


P9 


01:47:02 


01:47:12 


61±25 


105±18 


18.2 


P28 


01:48:37 


01:48:47 


42±42 


106±37 


9.9 


P30 


01:48:47 


01:48:57 


90±36 


54±11 


11.8 



Errors are given at 1 a c.l. for one parameter of inter- 
est. 

our first peak analysis, where we do not detect 
any polarized signal. By comparing P6 and PS, 
it seems that the polarization signal in this time 
interval results from the superposition of signals 
with different polarization angles, which could 
give rise to a null signal over a longer time in- 
terval. The shorter duration of the second pulse, 
where a polarization signal (P30) is clearly dom- 
inant over the others, could, on the other hand, 
explain the detection of a polarized signal over 
the whole second peak. Similarly, the 12 s peak 
emission may dominate the SPI data, while in the 
IBIS data this time period is heavily affected by 
telemetry losses, making the flux ratio between 
the brightest part (PS) and the rest of the GRB 
first peak smaller. In other words, the different 
polarization phases of the rest of the peak have a 
larger weight in our average measure. 

The expected level of polarization of the prompt 
7-ray emission in GRBs has been estimated by 
several authors for different models, or variations 
within them. In most cases, the observed 7-ray 
emission is due to the synchrotron radiation from 
relativistic electrons in the fast cooling regime. 
Their time-averaged distribution is a broken power 
law, 71(7) oc 7~^ with p' — p + 1 above Pm and 
p' — 2 below, where Pm is the minimum Lorentz 
factor of the injected dis tribution of electrons, 
and p ~ 2 — 2.5 its slope (jSari. Piran fc Naravan 
19981 ). The intrinsic polarization of the syn- 
chrotron radiation, Ilsvn = {p' + l)/(p' + 7/3) 
( Ribickv Sz Lightman 19791 ) is then of the order of 
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Ilsyn = {p + 2)/(p + 10/3) ~ 75% above and 
Ilsyn = 9/12 ~ 70% below, where Vm, the peak of 
the spectrum in vF^, is the synchrotron frequency 
of electrons at Fm. High polarization levels can 
also be reached if inverse Compton scatterings are 
the dominant radiative process. 

Our results show that (i) the polarization level 
in GRB04119A is varying on short time scales and 
can reach high values that correspond to a sizeable 
fraction of the intrinsic polarization Ilsyn. The 
polarization angle is varying as well, and (ii) the 
time-averaged value over longer intervals shows 
reduced polarization. We discuss these results in 
the context of different scenarios : 
(1) synchrotron emission from shock-accelerated 
electrons in a relativistic jet with an ordered mag- 
netic field contained in the plane perpendicular 
to the jet velocity. This geometry is favored if 
the field is carried by the outflow from the cen- 
tral source, as the poloidal component decreases 
much faster with radius than the toroidal one. 
The polarization level at the peak of a given pulse 



can be as high as n/IIsyn ~ 0.8, i.e. 11 ~ 60%, 
leading to a maximum time-averaged polarization 
in long intervals of ll/IIm ax ^ 0.6, i.e. 11 ~ 45% 
in this case (Granot 2003 : iGranot fc Konigl l2003t 
Nakar. Piran & WaxmanI l2003h . The main re- 
quirement is to have a uniform magnetic field 
in space, i.e. with a coherent scale R9b with 
}^ 1/r. The fact that the polarization level 
and angle vary during the burst indicates on the 
other hand that the field is not necessarily uni- 
form in time, as it would be required to explain 
a high level of th e timc-intcgratcd polarization 
([Nakar. Piran fc W axmau 2003). A magnetic field 
anchored in the central engine and c arried by the 
outflo w to large distance (see e.g. ISpruit et al 



20011 ) is therefore not the only possibility. A 
magnetic field generated at the shock could also 
work, and even favor variability, if there is a pro- 
cess capable to increase the field coherence scale 
(the field is most probably initially generated on 
small, skin-depth, scales). The existence of such 
a process is unclear in our present knowledge of 
the micro-physics in mi ldly relativistic shocks (see 
e.g. iKeshet et al. I l2008l) . Note that the condition 
^ 1/r is really necessary only in the pulses 
with the highest level of polarization. If Ob is 
smaller, so that a number N (F^b)^^ of mutu- 



ally incoherent patches are present is the visible 
region, the level of polarization will decrease, but 
the variability (both of the polarization level and 
angle) will increase (lGranotll2003h . If the radi- 
ating ele ctrons are accel erated in internal shocks 
(iRees fc JVIcsz aros 1 993: Kobavashi. Piran fc Sari 
19971 : iDaigne fc MochkovitchI Il998l ). the Lorentz 
factor is necessarily varying in the outflow, which 
can be an additional source of variability for the 
polarization. If 6b and l/F are close, the number 
of coherent patches in the visible region could vary 
from a pulse to another. We therefore conclude 
that any scenario where the observed gamma-rays 
are produced by synchrotron radiation from elec- 
trons in a relativistic jet with an ordered magnetic 
field in the plane perpendicular to the jet velocity 
seems fully consistent with our observations, as 
long as the coherence scale of the field is larger 
than l/F in most of the emitting regions. A poten- 
tial difficulty remains: an additional random com- 
ponent of the magnetic field is probably necessary 
to allow for particle acceleration in shocks. This 
component would reduce the level of polarization 
by a factor that is however difficult to estimate, 
as the int ensity of the the former is not well con- 
strain ed ( GranotI 20031 : Nakar. Piran fc Waxman 
2003h : 



(2) synchrotron emission from a purely electro- 
magnetic outflow. The estimated level of polar- 
i zation is comparable with the previous scenario 



(|Lvutikov et al. I I2OO3I) . In addition, a magnetic 



field with a large coherence scale is naturally ex- 
pected in such a purely electromagnetic outfiow. 
One potential difficulty is, however, related to the 
mechanism responsible for the energy dissipation. 
It this scenario, the energy has to be extracted 
from the magnetic field before being radiated. 
Therefore magnetic dissipation must occur in 
the emitting region, changing the field geometry, 
which becomes probably much less ordered, reduc- 
i ng the final level of po l arization by a large factor 
(iLvutikov et al. 1 I2OO3I: iNakar. Piran fc WaxmanI 
20031 ) ■ This effect is however difficult to estimate, 
as the details of the physical processes that could 
lead to magnetic dissipation in such an outflow 
are still far from being understood; 

(3) synchrotron emission from shock-accelerated 
electrons in a relativistic jet with a random fleld 
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generated at the shock and contained in the 
plane perpendicular to the jet velocity. A high 
level of polarization can be obtained even with 
a random magnetic field if the jet is observed 
from just outside its edge (jGhisellini fc Lazzati 
19991 : Waxman 2003 ). The polarization at the 



peak of a given pulse can reach II/II. 



syn 



0. 



I.e. 



n 



n 



60% resulting in a time-integrated 
0.5 - 0.6, 



40 - 45% (Granod 2003: 


Granot & Konig] 


Nakar. Piran & Waxman 


200d). However 



these high values are obtained if the jet is seen 
with 9ohs — + 1/r, where 9j is the opening angle 
of the jet and 6'obs the angle between the line-of- 
sight and the jet axis. Such viewing conditions 
are rare, except if 6j ^ l/F. Variability of the 
polarization level is expected if the Lorentz factor 
is varying in the outflow, as for instance in the 
internal shock model. Different pulses in the light 
curve correspond to viewing angles 6'obs = Oj+k/T, 
where k is larger for emitting regions with a larger 
Lorentz factor. The highest polarization is ob- 
tained in pulses with fc ~ 1 whereas the pulse 
flux decreases with fc for fc > 0. The highest po- 
larization should therefore not be found in the 
brightest pulses, which is difficult to test as all 
pulses do not have necessarily the same intrinsic 
luminosity. In addition, when several pulses super- 
impose, the measured polarization level, which is 
flux -weighted, could be red uced by a sizeable fac- 
tor ( Granot fc Konig|2003l ). Finally, the observed 



polarization can also b e reduced if the jet edges 



are n ot sharp enough (jNakar. Piran fc Waxman 
l2003[ ). We conclude that this model cannot be 
rejected as it can in principle reach high levels of 
polarization, but the conditions necessary for it 
seem more difficult to be achieved than in model 



(1). As shown bv I Granot fc Konig] (2003), similar 
conditions as for scenario (3) would be required 
for a scenario where the field is ordered but par- 
allel to the jet, leading to the same conclusions. 

(4) inverse Compton emission from relativistic 
electrons in a jet propagating within a photon 
field (" Compton drag" model) . The level of polar- 
ization in this scenario can be even higher than for 
the synchrotron radiation and reach 60 — 100%, 
but only under the condition tha t th e jet is nar- 
row with T9j < 5 (|Lazzati et al.ll2003 ). The max- 
imum level of polarization is again obtained for 



^obs — + 1/r- These viewing conditions are 
very similar to those of models (3) and are dis- 
cussed below in the case of GRB 041219A. Again, 
the polarization is reduced if the edges of the jet 
are not sharp enough. Variability of the Lorentz 
factor will again result in a varying polarization, 
with the same difficulties regarding the final level 
of polarization than in model (3). However, vari- 
ations of the Lorentz factor could possibly be less 
large in this scenario as part of the variability of 
the light curve can be related to the inhomogene- 
ity of the ambient photon field. 



Combining the Yonetoku relation (jYonetoku et al 



20041 ) to estimate the redshift and the isotropic en- 
er gy, with the standard energy reservoir d erived 
bv lFrafl et all (|2000h , iMcGlvnn et"all (|2007t l have 



estimated that 9j ~ 2-3 degrees in GRB 041219A, 
which leads to T0j - (3.5 - 5) (F/lOO). Such an 
estimate, based on very debated relations, indi- 
cates that the peculiar viewing conditions nec- 
essary in the geometric models are difficult in 
the case of GRB 041219A but cannot be fully 
excluded. Models (3-4) can therefore not be re- 
jected on the basis of these conditions only. As 
explained above, it is however not clear if a high 
polarization level can be maintained together with 
a high time variability. The status of model (2) 
is not clear as the details of the magnetic dissipa- 
tion process are crucial and not well understood. 
Finally, our results favor model (1) as it seems 
to have the best capabilities to maintain a high 
polarization level with some variability. An ad- 
ditional random component of the magnetic field 
is however necessary for particle acceleration at 
the shock and its negative impact on the level of 
polarization should be carefully estimated. A final 
answer to distinguish between intrinsic and geo- 
metric models could be obtained by accumulating 
more observations. Indeed, models (1-2) predict a 
polarized emission for all bursts, whereas models 
(3-4) would predict that only a small fraction of 
GRBs are highly polarized. This idea has been 
tested rec ently by M ontecarlo simulations per- 
formed bv iToma et a l. (2008), who conclude that 
if more than 30 % of bursts are polarized, geo- 
metric models can be ruled out. This shows the 
importance of polarimetric measurements for the 
understanding of intrinsic properties of GRBs, but 
the current instrumentation is statistically limited 
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and can provide measurements just for the bright- 
est events. To significantly increase the sample 
of GRBs with measured polarization, future GRB 
mission with wide field of views and polarimetric 
capabilities are needed. 
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